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The  spectrum  of  early  time  pressure  perturbations,  due  to  optically  smoothed  in¬ 
duced  spatial  incoherence  (ISI)  laser  imprint,  is  computed  for  a  planar  target  using 
a  forced,  dissipative  model.  The  time-dependent  ISI  laser  deposition  is  computed 
using  a  time-dependent  electromagnetic  full  wave  Maxwell  code.  It  is  found  that 
the  pressure  spectrum  evolves  into  a  power  law  in  which  spectral  power  is  trans¬ 
ferred  from  large  to  smaller  scales  through  a  nonlinear  cascade  process.  The  model 
results  are  compared  with  experimental  observations. 
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I.  INTRODUCTION 

A  limitation  on  high-gain  direct  drive  inertial  confinement  fusion  is  pellet  degra¬ 
dation  and  disruption  due  to  the  Rayleigh-Taylor  and  Richtmyer-Meshkov  instability1 
An  important  issue  is  the  nature  of  the  seeding  of  these  instabilities.  This  seed¬ 
ing  can  be  generated  by  pellet  surface  imperfections  and  laser  nonuniformities,  i.e., 
laser  imprint. 

Much  work,  both  experimental2-9  and  theoretical10-18,  has  been  performed  deal¬ 
ing  with  laser  imprint.  We  define  the  imprint  phase  as  the  period  before  the  return 
of  the  first  rarefaction  wave  to  the  ablation  surface  after  shock  breakout.  It  has 
been  demonstrated  that  the  size  of  the  conduction  zone  between  the  critical  surface 
and  the  ablation  surface  is  important  for  smoothing  laser  nonuniformities.  In  the 
start  up  phase  at  early  times  before  there  is  a  large  separation  between  the  critical 
and  ablation  surface  several  studies  have  shown  that  the  larger  scale  modes  seem 
to  dominate  the  smaller  scale  modes6,8.  To  minimize  the  effects  of  imprint,  much 
work  has  been  directed  at  improving  the  spatial  uniformity  of  the  laser  irradiance, 
e.g.,  using  optical  smoothing  methods  such  as  ISI19.  Some  nonuniformity  remains 
even  after  smoothing  takes  place. 

In  the  time-dependent  start-up  phase  the  pressure  profile  will  approximately 
match  the  driving  laser  intensity  profile.  For  the  seeding  problem,  it  is  important 
to  develop  a  model  for  the  spectrum  of  pressure  perturbations  over  a  wide  range 
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of  spatial  scales  at  early  times.  Nonlinear  effects  may  be  important  in  the  estab¬ 
lishment  of  the  early  time  pressure  spectrum.  In  this  study  we  develop  a  nonlinear 
model  of  the  pressure  perturbation  spectrum  due  to  optically  smoothed  induced 
spatial  incoherence  (ISI)  laser  imprint.  We  find  that  in  the  pressure  spectrum  a 
nonlinear  cascade  develops  where  energy  is  transferred  from  the  large  scales  down 
to  small  scales. 

II.  MODEL 

We  consider  a  planar  fully  ionized  plasma  with  the  initial  density  and  tem¬ 
perature  profiles  shown  in  Fig.  1.  The  laser  is  incident  from  the  left.  Using  a 
hydrodynamic  fluid  model  the  interaction  of  a  high  intensity  laser  with  a  planar 
target  can  be  described  using  a  complete  set  of  fluid  equations  which  are  coupled 
to  Maxwell’s  equations.  It  is  not  possible  to  solve  this  coupled  model  on  both  laser 
waveperiod  femtosecond  and  hydrodynamic  nanosecond  time  scales  simultaneously. 
Instead,  we  adopt  a  forced,  dissipative  model  to  simulate  the  hydrodynamic  evo¬ 
lution.  The  forced,  dissipative  approach  has  been  widely  used  in  hydrodynamic 
turbulence  studies.  This  approach  allows  the  study  of  aspects  of  the  fundamental 
question  of  coupling  among  a  range  of  spatial  scales  and  the  spectrum  of  pres¬ 
sure  perturbations.  We  adopt  the  following  equations  for  density,  momentum,  and 
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where  p  =  pT  =  pe+pi ,  Wji  =  l/2(dVi/dxj  +  dVj/dxi),  A  =  V- V,  q  =  -kVT,  rj  is 
the  viscosity,  k  is  the  Spitzer  thermal  conductivity,  5ij  is  the  Kronecker  delta,  and  7 
is  the  specific  heat  ratio.  Here  J  -E  gives  the  laser  energy  deposition  due  to  collisional 
and  resonance  absorption.  Using  a  collisional  fluid  model,  J  •  E  ~  oyE2  ~  arI  with 
a  =  iuipe/(4:iT(u}  +  ivei))  where  upe  is  the  plasma  frequency,  u>  is  the  laser  frequency, 
uei  is  the  collision  frequency,  and  I  is  the  laser  intensity  By  scaling  the  pressure, 
velocity,  length,  and  temperature  by  po,  bo,  L0,  and  T0,  Eq.  (l)-(3)  give  three 
nondimensional  parameters,  i.e.,  a  Reynolds  number  Re  =  poVoL0/r],  Vo/Cs,  and 
tj/k  where  Cs  is  the  sound  speed. 

We  are  interested  in  developing  a  model  for  the  early  time  start-up  phase 
pressure  perturbation  spectra  from  lasers  smoothed  by  ISI  methods.  In  a  two- 
dimensional  plane  perpendicular  to  the  direction  of  laser  propagation,  with  / (x,  y )  = 
E*  •  E.  an  expression  for  the  spectrum  of  laser  intensity  using  ISI  smoothing  can  be 
written16: 

I  ( k )  =  ^7r  —  2sin~1k  —  2k  (l  —  k2^j  (4) 


2  Nz 


where  k  =  {k2x  +  /c2)1/2,  Nz  =  kmaxDspot/2n,  krnax  =  27t//A0,  with  Dspot  the  laser 
spot  size,  f  is  the  f  number,  Aq  is  the  laser  wavelength,  and  Iavg  is  the  average 
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laser  intensity.  We  note  that  the  ISI  smoothed  laser  intensity  deposit  energy  over  a 
range  of  spatial  scale  sizes.  As  a  result,  from  Eq.  (3)  to  lowest  order,  the  pressure 
perturbations  will  also  be  driven  over  a  range  of  spatial  scale  sizes.  Nonlinear  effects 
may  couple  these  spatial  scales.  Knowing  the  intensity  spectrum  I(x,y,z)  the  full 
spectral  content  of  J  •  E  from  an  ISI  smoothed  laser  needs  to  be  calculated  in  order 
to  solve  Eq.  (3).  We  have  computed  the  spectrum  of  J  •  E  for  ISI  laser  intensity 
of  10  -  50  TW/cm2  and  over  a  range  of  scales  using  a  time-dependent  Maxwell 
equation  full  wave  code: 


V  x  E 
V  x  B  = 


ldB 

c  dt 

4tt  10E 
c  ^  c  dt 


J  =  en(Vi-Ve) 


(5) 

(6) 
(7) 


where  Ve  and  V,  are  found  from  the  electron  and  ion  momentum  equations.  Fig. 
2  shows  the  k-spectrum  of  J  •  E  using  the  full  wave  Maxwell  code.  The  laser  was 
taken  to  be  normally  incident  on  the  plasma  slab  configuration  shown  in  Fig.  1. 

We  solve  Eq.  (l)-(3)  numerically  using  a  pseudospectral  method.  All  transverse 
spatial  derivatives  are  solved  in  Fourier  k-spaee.  The  Eve  independent  variables  are 
P,P,Vx,Vy,Vz.  The  time  step  algorithm  is  performed  by  a  Runge-Kutta  scheme.  We 
use  128  x  128  x  64  grid  points  in  the  x,y,z  directions,  respectively.  The  boundary 
conditions  are  taken  to  be  periodic  in  the  transverse  (x,y)  directions  and  outflow 
in  the  z-direction.  The  simulation  volume  has  dimensions  of  60  /jrn,  60  /irn,  and  10 
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/jrn  in  the  x,y,z  directions.  As  initial  conditions,  our  approach  is  to  force  Eq.  (3) 
using  spatial  modes  taken  from  the  spectrum  of  J  •  E  as  computed  near  the  critical 
surface  using  the  full- wave  Maxwell  code  of  Eq.  (5)-(7).  We  then  follow  the  forced 
nonlinear  hydrodynamic  evolution  and  study  the  nonlinear  coupling  between  the 
forced  modes  and  the  smaller  scale  modes.  Aliasing  errors  are  treated  by  truncating 
the  Fourier  values  outside  a  shell  of  hxed  wavenumber  kmax  =  21^2N/3  where  N  is 
the  number  of  grid  points  in  each  direction.  We  take  the  following  form  for  the 
forcing  function: 

S  —  J  ■  E  =  So  (A  (t)  sink  ■  x  +  B  (t)  cosk  •  x)  (8) 

where  A  and  B  are  random  variables.  Each  element  A  and  B  is  assumed  to  be 
statistically  independent.  The  absolute  magnitudes  for  So,  A,  and  B  are  derived 
from  the  J  •  E  spectrum  computed  from  Eq.  (5)- (7).  In  this  study  we  force  only 
the  spatial  modes  k  with  27t/50 pm  <  |k|  <  2n/ 40  pm.  This  will  allow  the  study 
of  the  cascading  effects  down  to  smaller  scale  modes  to  compare  with  experimental 
observations. 

III.  RESULTS 

By  defining  the  total  energy  by  Et  =  Ej  +  Ek  with  Ej  =  p/(y  —  1)  and 
Er  =  (p/ 2)V2,  we  have  found  that  the  total  and  internal  energy  increase  in  time 
due  to  the  ablation  pressure  applied  by  the  J  •  E  laser  intensity  profile  but  the 
kinetic  energy  stops  increasing  and  oscillates  about  a  quasiequilibrium  value.  We 
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have  also  found  that  Ek  is  much  smaller  than  Ej  for  the  cases  studied. 

Fig.  3  shows  the  time  development  of  the  maximum  and  spatially  averaged 
pressure.  After  an  initial  transient  period  the  average  pressure  slowly  increases 
linearly  in  time.  Fig.  4  shows  a  plot  of  the  pressure  contours  and  illustrates  the 
transverse  (x,y)  structure.  Fig.  5  shows  the  time-averaged  pressure  spectrum.  We 
have  defined  the  pressure  spectrum: 


1  /  P_ 

2  \Po , 


=  j  P(k)  dk 


(9) 


with  p  =  p  —  p0.  We  have  found  that,  for  the  cases  studied,  the  presssure  spec¬ 
trum,  from  the  simulation  model,  can  be  fit  with  a  power  law  of  the  form  P{k)  ~ 
k~lexp(—km)  for  the  approximate  scale  size  regime  of  60  /jrn  -  3  /jrn. 

An  approximate  analytical  expression  for  the  pressure  spectrum  can  be  obtained 
as  follows.  In  Eq.  (3)  we  let  Vj  =  Vq  +  Vj,  S  =  So  +  Sj,  and  write: 


dt  2 


TTy  dp  dp  dp  Tj 


-pVi 
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—  K 


dxj  dxj 


+  pS 
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where  /  denotes  time  average.  We  then  take: 


l(^)2  =  IE^dk 


(ii) 


We  find  that  the  pressure  slowly  increases  in  time  reaching  a  quasiequilibrium.  We 
take  as  an  approximation  d/dt  ~  0  and  balance  the  third  and  fourth  terms  on  the 


right  hand  side  and  find: 


P  (k)  ~  k  Pexp(—Akq ) 


(12) 


Here  we  have  let  E(k)  —  Eok~n  with  p=2+(l-n)/2  and  q=(n+l)/2.  From  the 
cases  studied  we  find  n  ~  1.5  —  2,  p  ~  1.5  —  1.8,  and  g  ~  1.3  —  1.5.  We  have  found 
that  the  pressure  spectrum  is  approximtely  isotropic.  This  can  be  seen  as  follows. 
Eliminating  the  velocity  from  Eq.  (2)  and  inserting  into  the  pressure  equation  Eq. 
(3)  we  can  write  in  k,  tv  space,  after  ignoring  source  and  sink  terms,  the  lowest  order 
nonlinear  equation: 


9pk,u>  _  \  '  rp  ~ 

/  j  ^  k,k;  k7  ,u>— ui' 

k/cj/ 


(13) 


with 

1  i  (k  •  k'  +  (7  —  1)  /c'2) 

Tk’kW  =  y0  u  -  i^k’2 

and  po  is  an  average  density.  We  note  that  in  the  expression  for  the  nonlinear 
matrix  element  T,  which  describes  the  nonlinear  coupling  to  lowest  order,  for  7  >  1 
the  second  term  is  larger.  The  second  term  is  function  of  k '  =  |k  |  suggesting  an 
approximately  isotropic  nonlinear  coupling  to  lowest  order. 

We  have  found  that  the  spectrum  of  pressure  perturbations  develops  into  a  power 
law  where  large  scales  are  dominant  over  smaller  scales.  Previous  experimental 
studies  in  the  start-up  phase6’8,  using  similar  laser  intensities  and  scale  size  regimes, 
have  also  found  a  power  law  dependence  in  k-space. 
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IV.  SUMMARY 

In  summary,  we  have  developed  a  nonlinear  model  for  the  spectrum  of  pressure 
perturbations  due  to  optically  smoothed  lasers.  This  has  been  accomplished  using  a 
forced,  dissipative  model.  The  time  dependent  laser  deposition  has  been  computed 
for  ISI  using  a  full  wave  electromagnetic  Maxwell  code.  We  have  found  that  the  early 
time  pressure  spectra  develops  into  a  power  law  dependence  where  large  scales  are 
dominant  over  smaller  scales.  We  have  found  that  the  model  results  are  consistent 
with  previous  experimental  observations.  In  the  future  we  will  study  the  pressure 
spectra  from  other  spatial  scale  size  forcings. 

This  work  was  supported  by  the  Department  of  Energy.  We  thank  Steve  Oben- 
schain,  A.  Vclikovich,  D.  Fyfe,  J.  Gardner,  S.  Zalesak,  D.  Colombant,  and  W. 


Manheimer  for  useful  discussions. 
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Figure  Captions 

Fig.  1  Plot  of  initial  density  N  and  temperature  T  configuration. 

Fig.  2  Plot  of  angle-averaged  k-spectrum  of  J  •  E.  Normalizations  of  Sq  =  1023  erg 
cm-3  sec-1  and  ko  =  27t/60 fim. 

Fig.  3  Plot  of  time  history  of  pressure.  The  pressure  has  been  normalized  by 
P0  =  4Mbar  and  to  =0.15  psec. 

Fig.  4  Plot  of  normalized  transverse  pressure  profile.  The  pressure  here  has  been 
normalized  by  5  Mbar.  Eight  evenly  spaced  contours  have  been  plotted. 

Fig.  5  Plot  of  pressure  power  spectra.  Here  Pq  =  25Mbar2  and  ko  =  2n/60 fim. 


Fig.  2 


CL|  Q. 


